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We calculate the relative conserved currents, superpotentials and conserved quantities between two 
O \ homogeneous and isotropic universes. In particular we prove that their relative "energy" (defined 

■ as the conserved quantity associated to cosmic time coordinate translations for a comoving ob- 

' server) is vanishing and so are the other conserved quantities related to a Lie subalgebra of vector 

fields isomorphic to the Poincare algebra. These quantities are also conserved in time. We also 
find a relative conserved quantity for such a kind of solutions which is conserved in time though 
non-vanishing. This example provides at least two insights in the theory of conserved quantities 
in General Relativity. First, the contribution of the cosmological matter fiuid to the conserved 
O quantities is carefully studied and proved to be vanishing. Second, we explicitly show that our 

superpotential (that happens to coincide with the so-called KBL potential although it is generated 
»' I ' differently) provides strong conservation laws under much weaker hypotheses than the ones usually 

{^X), required. In particular, the symmetry generator is not needed to be Killing (nor Killing of the 

background, nor asymptotically Killing), the prescription is quasi-local and it works fine in a finite 
region too and no matching condition on the boundary is required. 
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I. INTRODUCTION 



The definition of energy is one of the most disputed topics in General Relativity. In the case of isolated systems, 
such as black holes, different prescriptions are available to define energy (as well as for the other conserved quantities) 
, and on the simplest known examples where most of them are applicable, they usually provide the same result [3] . In 
' these prescriptions at least one of the following requirements is always asked for: the symmetry generators are Killing 
. !^ \ vectors (at least of the background, or at least asymptotically) ; some matching between the dynamical metric and 
. the background; some special asymptotic behavior. In the case of cosmological solutions, as we shall see, all these 
$H ' techniques cannot be applied and there is no general agreement not only on how much energy is there (nor in the 
universe nor in its local domains), but we do not even have a clear definition of what the energy should be. There 
are References (see for example Refs. [4, 30, 34]) where the energy of some models of homogeneous and isotropic 
universes is defined according to different prescriptions (mainly relying on pseudotensors) and it turns out to be 
vanishing; however some criticism to these non-covariant methods can be found for example in Ref [8] and General 
Relativity principles are quite clear on this: either such quantities can be described in some covariant way or they 
are not fundamental to the description of the Universe. Recently three articles have been published with definitions 
leading to different results (see. Refs [2, 21, 22]). 
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From the purely physical viewpoint the absence of a region far off from the sources makes very hard to give a physical 
interpretation to conserved quantities that however can be calculated on a purely mathematical ground under very 
general hypothesis once one clearly fixes the definitions. 

In recent years many authors ([1, 23, 24, 25, 31]) provided prescriptions to compute the conserved currents in 
General Relativity based on a suitable "boundary correction" to the Komar superpotential relative to a reference 
configuration. 

To this research effort two of the present authors gave a number of contributions (see for example Refs.[ll, 13, 15, 
16, 17]) mainly proposing a pure variational route to the definition of the relative conserved currents between two 
solutions of the field equations interpreted as the amount of conserved current needed to pass from one solution to 
the other. The mathematical framework introduced to this aim goes well beyond General Relativity and extends in 
fact to any gauge- natural Lagrangian field theory (see for example Ref. [14] and Ref. [9]). 

The aim of this paper is to apply and adapt the previously known techniques to the calculation of the relative 
conserved quantities in the case of two isotropic cosmological solutions with different densities, where some technical 
requirements commonly found in the literature are not fulfilled and where the contribution of matter has to be carefully 
taken into account. 

The physical motivation behind this choice is to understand how much a change of the scale factor (and thence of 
the matter amount of the Universe) affects the conserved currents and in particular the gravitational energy. 

We find that all the relative covariantly conserved quantities related to a Lie subalgebra of the Lie algebra of 
vector fields, isomorphic to the Lie algebra of the Poincarc group, are vanishing and in particular this holds true for 
the relative energy that we define as the conserved quantity associated to cosmic time translations for a comoving 
observer. Let us remark that these observers are not chosen a priori (this would not be allowed by the covariance 
principles of General Relativity) but they are a posteriori selected by this family of cosmological solutions as preferred 
cosmological observers who provide a better interpretation of the results. In principle our prescription applies to any 
observer. According to our interpretation we can say that "injecting" matter into the Universe do not change its 
conserved quantities. 



II. THE VARIATIONAL FRAMEWORK FOR RELATIVISTIC HYDRODYNAMICS 



Let us review here a variational formulation of relativistic hydrodynamics that was originally presented in detail 
in Ref. [5]. It is in practice equivalent to the ones of Refs. [26, 29, 32] but it shows some formal and computational 
advantages. The aim of reviewing it here is to prove (and not just assume as it is usually done) that the matter 
Lagrangian does not contribute to the conserved quantities related to the diffeomorphism invariance of the theory. 

Our kinematical description is given trough the following fields: gravitation is represented by a Lorentzian metric 
with components g^i^, while the fluid degrees of freedom are encoded in a nowhere vanishing vector density J = J^ds^ 
(where ds^ is the natural basis of (to — l)-forms over an TO-dimantional manifold. The form J is such that the unit 
timelike vector = J'^ /\J\ is tangent to the flow lines and that the volume form \/[gjp = | J| represents the matter 
density. Here g denotes the determinant of the metric tensor. 

In order to implement conservation of matter we only allow for closed J (i.e. d^J^ = = J^), so that for any 
closed 3-surface S in spacetime the flow of J through S vanishes. 

Starting from the Lagrangian 

L = Lh + Lf (1) 

with 

Lnigap, dxgap, dx^gafj) = {R - 2A) ds 

LFiJ"", gap) = -^M [P (1 + e(p))] ds = -^\^i{p) ds 

where R is the scalar curvature of the metric. Lambda a constant and e{p) a suitable function (here generic) of the 

scalar p = \l ^'^"\g\ '^ '~ ^ physically interpreted as the internal energy of the fluid. Pressure is defined as P = p^f|- 

We want to minimize the action functional preserving the constraint 9q, J" = 0. To do so we take arbitrary variations 
for the metric, while J is varied according to the rule 5J^^ = £xJ^ = d^J^X'' - J'^d^X'^ + J'^d^X^. Intuitively, 
it amounts to drag the flow lines of the fluid in the spacetime along an arbitrary diffeomorphism (of course fixing 
boundary values). This euristic prescription can be formally justified on rigorous bases; see Ref. [5]. 
The Euler-Lagrange equations arising from this "constrained" prescription are the following 
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V,, = 



(constraint) 



where it^ = J'^/lJ] and H^^ is the fluid stress tensor 



(2) 



Pdt^u + + P)U'U 



DifFeomorphism covariance of the Lagrangian (1) provides the following Nothcr current associated to any spacetime 
vector field ^ = ^'^di, 



£{L, i) = i) ds^ = £h{Lh,S,T dso. - V\9\ (s"'' ^a-D ds^ 



where we set 



2k 



(3) 



(4) 



for the Nother current obtained from the Hilbert Lagrangian Lh- The current £{L,^) uniquely splits according to 
the general theory on gauge-natural field theories (see [14], chapter 6) into = £{L,^) + Div U{L,S^) where 

the reduced current £{L,^) vanishes on shell (i.e. when computed along a solution of the field equations) while the 
superpotential U{L,^) has the usual Komar expression 



zk 



(5) 



The fluid Lagrangian does not contribute to the superpotential since its contribution to the Nother current is of order 
zero in ^ (see again Ref. [5] and the general theory in Ref. [14] for more details) . 

III. RELATIVE CONSERVATION LAWS FOR HOMOGENEOUS AND ISOTROPIC UNIVERSES AND 

RELATIVISTIC HYDRODYNAMICS 



A homogeneous and isotropic cosmological model is a particular class of solutions of the system described by the 
Lagrangian (1) where the couple (.g, J) is such that in a system of synchronous comoving coordinates (t,r, 0, (/>) it 
looks locally like (see for example [28]) 



1 - A>2 



and 



J 



where dso = dt-ids and K G {—1,0, 1} is the normalized constant scalar curvature of the spatial leafs at constant t, 
while the scale factor a{t) and the density p are subjected to equations (2) that in the chosen coordinates read as 



Sct^ + 3fc - Kfia^ - Aa^ = 



(6) 



being in this case the third equation in (2) equivalent to the first one. 

Quasi-local conserved quantities in a closed region D of spacetime for a given solution of the field equations and 
relative to a Lagrangian symmetry can be defined as the integral of the superpotential along the boundary of D. This 
definition however is affected by a number of ambiguities that have been thoroughly discussed in the literature (see 
for example Ref. [14] Chapter 6, Section 5 and Ref. [11]). 

A better defined concept is that of relative conserved currents and quantities. They arise in any gauge- natural 
Lagrangian field theory (of any order) from the following prescription. Let be j/* be the fields, their derivatives (first 
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order for simplicity, but higher order are allowed too) and L{x^ , , y^^) a gauge-natural Lagrangian (any Lagrangian 
that is covariant under spacetime difFeomorphisms or gauge transformations belongs to this class, cf. Ref. [14]) with 
first variation formula 5y-L = E{L, Y) + Div F(L, Y) where E{L, Y) = Q are the field equations and F{L, Y) is the 
so-called Poincare-Cartan morphism. Let us consider a one parameter family {2/^('^^)}se[o,i] of solutions of the field 
equations such that j/g = is interpreted as a reference background and y\ = j/Ms a solution representing the physical 
fields. Our prescription amounts to correct the superpotential in order to fix the zero level for the conserved quantities 
in such a way that they become vanishing when the physical fields coincide with the given reference. To this aim, 
however, covariant- ADM formalism (and also other derivations, sec Refs.[13, 18, 19, 20] ) suggests to consider a new 
relative superpotential U'^{^,ys,y) built from the difference of the values of the old superpotential along yl and y*, 
corrected in such a way that its variation along the infinitesimal generator X of the 1-parameter family {j/s(2;^)}se[o,i] 
gives 

SxU'i^, ys, y) - 5xU{^, y) + ^F{L, X). (7) 
In vacuum General Relativity a superpotential the variation of which satisfies (7) is known to be the following one 

Whi^. g, -g) = Vl-e'^l - /^Vl-^'^l + .Mr'-^t^'^^ds^P- (8) 

where g is interpreted as the physical metric while g is the reference configuration (sometimes called a dynamical 
background as it is required to be a solution of Einstein equations) and = F^^ — f^^^F^j^ — f — (^^^f ^j^. Barred 
quantities are defined with respect to the background metric. 

This relative superpotential, as claimed, is the difference between the two copies of the Komar superpotential (5) 
evaluated on g and g plus the covariant- ADM correction and it vanishes when g = g. This superpotential happens to 
coincide with the so called KBL superpotential (named after the authors of Ref. [25]) that, however, was derived in 
a different ad-hoc way and under stricter hypothesis rather than being shown to be the specific case of a much wide 
concept. 

We define now 

Q{^,D,g,g)= [ U'{tg,g)- (9) 

JdD 

as the relative conserved quantities between the two (homotopical) solutions g and g in the domain D and relatively 
to the Lagrangian symmetry generator ^. In other words Q may be interpreted as the cost in terms of the conserved 
quantity needed to pass from the solution g to the solution g (for any fixed region and for any vector field) . 

We point out that besides many side conditions are usually required in the literature in order to give a precise 
physical interpretation to the conserved quantities (^ is often asked to be a Killing vector either of g or g, D is 
the spatial infinity, g and g are usually required to match on dD), since the Hilbert Lagrangian is invariant under 
diffeomorphisms, the relative superpotential (8) is a strongly conserved current associated to any vector field ^ without 
any restriction on ^ itself, on D, and also without boundary conditions on g and g. In Appendix A we review a further 
theoretical distinction between general conserved quantities wit respect to those that are also conserved in the time 
induced by an ADM foliation of spacetime. 

In presence of a cosmological fluid we already noticed that the matter contribution to the Komar superpotential 
is vanishing. Let us however explore the need for some covariant-ADM-like corrective boundary term such as that 
appearing in the relative superpotential (8) when one introduces as a reference background another FRW universe 
with a different amount of matter (thus different scale factor), but same topology for the spacial leaves. The choice 
of such a background is motivated by the idea of understanding how much a change of the scale factor (and thence of 
the matter amount of the Universe) affects the conserved currents and in particular the gravitational energy. 

Such a corrective term needs to satisfy condition (7) for the fluid Lagrangian that reads 

dxU^i^,gsrg) = SxiiFiLF,X) 

with X the generator X = {Xg,Xj) of the 1-parameter family of solutions {gjf} , J^^-^) with (.gfiV , >/("■)) being the 
physical configuration (a FRW universe with or without cosmological constant) and {gl?J , J"q-^) the reference one 
(another FRW universe with a different amoimt of matter, but same topology) and where F{Lf,X) is the term that 
goes under divergence in the first variation of the fluid Lagrangian (the so-called Poincare-Cartan morphism) taken 
with respect to X. 
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In order to find such a boundary correction to the superpotential let us notice that in a homogeneous and isotropic 
cosmological model the unique parameter left free by the symmetry requirements is the initial condition for the second 
of equations (6) and thence a 1-parameter family (gs, Jg) of homogeneous and isotropic universes is allowed to depend 
on the parameter just through the expansion factors as{t). Let us denote from now on by a{t) and a{t) respectively 
the expansion factors of g and of g. Due to the first of equations (6), the J, cannot depend on as{t) and thence they 
are all the same Js ^ J for each s. The second component Xj of the generator of the family is henceforth vanishing. 
As the fluid Lagrangian does not involve derivatives of g, the candidate corrective term is necessarily vanishing. 

In the case of homogeneous and isotropic cosmological solution, thence, the superpotential associated to any in- 
finitesimal generator ^ of diffomorphisms has still the same expression as that of equation (8). 

No restriction on ^, however, is required for the conservation law V^f = to hold on shell, nor for the on-shell 
identity £^ = V^Z-/'"^! and thence nor for the the strong (off-shell) conservation of Vi^Z^/I^'^l. 

In relation with the theory proposed in Ref. [9] we remark that the superpotential for the isotropic cosmologies can 
be directly derived from the augmented Lagrangian 

Lhf.us = Luig) - Lnig) - d^i^l^" W^^) + Lp{J,g) - LpiJ^g) (10) 

This Lagrangian is invariant under difFeomorphisms. To any vector field ^ we can associate the relative Nother 
current froi(i_ffFaug, ^, (ff, J), (.9, J)) = SnFrddsa where 

^HFrcl = ^HFia, J) - £%f{9, J) + dp [Mr^'v^tS.^^) (11) 

and £fjp{g, J) and £%p{g, J) can be calculated according to formula (4). 

Superpotential (8) arises from the splitting of the current (11) according to the general theory (see again Ref.[14], 
chapter 6). 

IV. THE VANISHING OF THE RELATIVE ENERGY AND DISCUSSIONS ON THE OTHER 

CONSERVED QUANTITIES 

Even if the definition of energy is not an undisputed topic, we already said that in many cases it has been fruitfully 
defined as the flow of the conserved current associated to the generator of time translations through a surface S at 
t = const delimited by a 2-dimensional boundary 91]. Cosmological solutions carry a preferred notion of time to be 
used to this aim. 

Let us remark that if t is the cosmic time, dt is not a Killing vector for g nor for g, thence these metrics are 
not invariant for cosmic time translations. However dt as any other generator of difFeomorphisms is a Lagrangian 
symmetry for (1) and this allow to calculate the strongly conserved current and quantity related to it by Nother 
theorem. 

To call the above quantity energy is now partially a matter of interpretation (and thus questionable in absence of 
a clear physical definition for this kind of spacetime); however, at least from the mathematical viewpoint, this is a 
well defined covariantly conserved current associated to a clearly distinguished time and we cannot see any strong 
argument against its interpretation as a kind of energy. Moreover, this is a trivial extension of the same techniques 
which have proven to be effective in the cases where the physical interpretation allows for a clear identification of 
what should be understood as energy on a physical instance (cf. [6, 7, 10, 11, 12, 13, 27]). 

Let us thence compute the relative energy between two homotopic homogeneous and isotropic universes (thus two 
universes with the same fT, no matter if it is —1,0 or 1, but different scales a(i) and a(t)) as the conserved quantity 
Q(Laug7 c?t, S, (5, J), (g, J)) with respect to a hypersurface E. 

Straightforward computations (here performed with the aid of Maple tensor package) lead us to conclude that the 
Nother current E = £^ds^ (11) has the first component vanishing on-shell, while the other components are identically 
vanishing. The superpotential is identically vanishing (in agreement with general theorems on the uniqueness of the 
splitting) and in particular all of its three summands vanish independently (the relative correction plays no role). 
Whichever is the surface E we conclude that the relative energy contained inside 9E is vanishing. From the physical 
viewpoint this means that not only the total energy of the universe as a whole vanishes, but also locally (quasi-locally) 
in any 3-dimensional domain bounded by a 2-surface 9E there is no total energy, no matter how much matter is there 
contained. A possible physical interpretation of this fact is that the internal energy of the cosmological fluid increases 
in fact with the density and it fuels the expansion, while the energy of the gravitational fleld increases with the density 
too but it fuels contraction (being a bonding interaction it contributes with the minus sign) and the two contributions 
balance giving rise to a vanishing total energy. Let us remark that the energy (being independent on S) besides being 
a covariantly conserved quantity is also in this case conserved in time. 
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Let us now consider the vector fields Tx,Ty,Tz defined as follows 



= sin 9 cos ( 



Ty = sin 9 sin c 



d cosScoso!) d sine 



dr r do r sin 9 84 

d cos 6* sin d cos0 d 



dr r 89 r sin 9 dcj) 



d sin9 d 
= cos6'— 



dr r 89 

These vectors are the generators of translations along the following cartesian-like coordinates x = r sm{9) cos{(f>) , 
y = rsm{9) sin((/)), z = rcos{9). Let moreover dTitQ,ro be the 2-dimensional hypersurface defined by r = ro and t = tg. 
The conserved quantities 

Q{Laug, Tx, S(t(,^ro)7 J)j (ff: J)) = 
= Q{Li,VLg,Ty, S(to,ro); (ff; J)^ id, J)) = 

Q(Laug,Tz,S(to,ro), (5, J)i idi J)) = 

are vanishing (and thus also conserved in time according to the definition of Appendix A). Let us however notice that 
in these cases the superpotentials are no longer vanishing and a different choice of the surfaces on which to compute 
the integral may lead to a different result. 

For the vector fields = yTz — zTy, Ly = zT^ — xT^ and = xTy — yT^ that generate rotations of the 
hypersurphaces of homogeneity we also have 

Q{Laug, Lx, S(t„^ro)7 {9, J), {9, J)) = 

= Q(Laug,iy,S(t^^rg), {g, J), {g, J)) = 

= Q(Laug,i2,S(to,ro)> (9, J)i (.9, J)) = 

again without an identically vanishing superpotential. 

If we now consider the vector fields ~ xdt + tT^ , By — ydt — tTy and B^ — zdt — tT^ that generate boosts in the 
coordinates {x,y,z) defined above we find 

Q{L&us, S3;, E(tg,r^), (g, J), (5, J)) ^ 

= QiLaug, By, Sft^ ,,0)7 (9, J)^ (ff: J)) = 

= Q(iaug,S3,I](t„_r„),(g, J),(5, J)) = 

again without an identically vanishing superpotential. 

The subalgebra of the Lie algebra of vector fields generated by the ten vector fields 9t, T.^, Ty, T^, L^, Ly, L^, B^, 
By, Bz is isomorphic to the Lie algebra of the Poincare group. The covariantly conserved quantities associated to 
them inside a surface dllto,ro defined by r = t-q and t = to are also conserved in time and vanishing. 

A natural question then arises also in consideration of the analysis carried out in a very famous paper on vacuum 
fluctuations (see Ref. [33]): is it possible to find a non- vanishing covariantly conserved quantity inside dT,to^roj 
possibly also conserved in time? The answer is yes. If we compute the conserved quantity associated to the vector 
field Z = dt for an homogeneous and isotropic universe with negative spatial curvature we find 

(9, J), (9, J)) 



4Vl + ^2 

that is also conserved in time (according to the definition of Appendix A). 

Let us notice that homogeneous and isotropic universes arc quite different from the other cases where the relative 
energy has been computed, at least for two reasons, one technical, one interpretative. From the technical viewpont, 
in fact, the generator of time translations is not a Killing vector for the dynamical metric nor for the background, 
and Dirichlet boundary conditions are not satisfied. We stress however that these conditions are not necessary for the 
conservation of the energy: in our case, in fact, they are not fulfilled but the current (11) is still conserved on-shell and 
the integral of the relative superpotential (8) is strongly conserved. Let us moreover remark that for homogeneous and 
isotropic universes it is not possible to ask for a matching of the dynamical metric with the background on dD since 
such a requirement would imply g — g. Nevertheless the Nother current can be computed with two non matching 
g and g and it is conserved. The unfulfillment of Dirichlet boundary conditions does not affect the variation of the 
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supcrpotential (8) that again provides (7). The covariantly conserved quantities associated to the generators of the 
Poincare subalgebra are moreover conserved also in time (according to the definition of Appendix A). 

In the previously known standard examples the matching on the boundary, or the Killing property for the symmetry 
generator, were always explicitly asked for in order to ensure the vanishing of the integral (A2). This is, in our 
experience, the first example where both requirements are not satisfied still giving rise to time conservation (for all 
the symmetry generators considered). 

From the viewpoint of physical interpretation we have to emphasize that, as was already pointed out by [8], isotropic 
universes are not asymptotically flat, and there is no way of going " infinitely" far away from the sources to move a 
the test particle and check if our definition of "energy" is compliant with intuition. In any case our definition provides 
a covariantly conserved quantity associated to time translations that is also conserved in time and that certainly has 
a physical meaning. 

Let us also remark that even if our starting point in deriving the superpotential (8) was to consider a 1-parameter 
family of solutions {gs} connecting g and g, nothing prevent us from considering from the very beginning the La- 
grangian (10) with two solutions that are not homotopic (of course in this case it is no longer possible to study the 
variation of its superpotential, and the covariant ADM formalism is not any longer a valid motivation, but let us go on 
and se what comes out). In the case of two homogeneous ad isotropic cosmological solutions with different curvature 
(e.g., —1 for the dynamical solution, for the background) we get that the relative energy inside the surface t — const, 
r = const is 



and considering the limit for r —f oo we get an infinite relative energy between the two configurations. This result is 
not surprising since to switch from a universe with flat hypersurfaces of homogeneity to another one in which they 
have negative curvature, a change in the topology is needed and it seems reasonable that one needs to spend an 
infinite energy to perform such a task. We stress that we do not have a proof of the fact that energy diverges if and 
only if g and g are not homotopic though to the best of our knowledge the only cases in which the energy is known 
to diverge are those when g and g are not homotopic. 



We applied variational calculus to calculate the relative covariantly conserved quantities between two homogeneous 
and isotropic cosmological models. We found that if the two models share the same normalized spatial curvature their 
relative conserved quantities relative to an entire subalgebra of vector fields isomorphic to that of the Poincare group 
are vanishing. In particular their relative energy, defined as the conserved quantity associated to time translations for 
a comoving observer, is vanishing. We also found that it is possible to chose a vector field such as that the associated 
conserved quantity is both non vanishing and conserved in time. 

Now that we extended the prescription of Ref. [9] to constrained systems (and in particular to relativistic fluids), 
further investigations will be devoted to describe conserved quantities of other systems of the same kind such as 
Tolman-Bondi collapse solutions. 



In the text we illustrated what we mean when saying that a quantity is covariantly conserved. In the literature, 
however, the word conserved is often associated to quantities that do not change in time. This is a different property 
and it is not manifestly covariant. However let us explain in which sense some covariantly conserved quantities are 
also conserved in time. Let us consider a spacelike 3-dimensional hypersurface S and a timelike vector field Y . Let 
us assume that the flow of Y drags the surface E generating an ADM foliation Es where s is the time associated to 
it at least locally. One can consider the 4-dimentional region D spanned by the surfaces Es when the parameter runs 
from So to si. It is bounded by the three surfaces E^^, Esj^ and B = UgSEg. Let us then consider a solution y of 
field equations, a symmetry ^ of the Lagrangian L and the conserved current £(L,^, y). The covariantly conserved 
quantity Q{L, ^, E^g , y) obtained by integrating £{L, ^, y) along Es„ is also conserved in time s if such an integral does 
not depend on which of the surfaces E^ the integral is computed on. Being dD = E^,, U E^^ U B one has 




a^{t)r^ 
2a^{t) 



V. 



CONCLUSIONS 



Appendix A: CONSERVATION LAWS "IN TIME' 



JD JdD JSt-^ J^ta JB 



(Al) 



8 



and thus conservation in time is equivalent to the vanishing of the flow of the current along B: 

I £iL,^,y)^0. (A2) 
Jb 

Different sets of sufficient conditions can be given in order to ensure the conservation in time of a covariantly conserved 
quantity. We do not list them here and we refer the reader to Ref. [13]. 
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